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Abstract 

We completely characterize General Gauge Mediation (GGM) at the weak scale 
by solving all IR constraints over the full parameter space. This is made possi¬ 
ble through a combination of numerical and analytical methods, based on a set of 
algebraic relations among the IR soft masses derived from the GGM boundary con¬ 
ditions in the UV. We show how tensions between just a few constraints determine 
the boundaries of the parameter space: electroweak symmetry breaking (EWSB), 
the Higgs mass, slepton tachyons, and left-handed stop/sbottom tachyons. While 
these constraints allow the left-handed squarks to be arbitrarily light, they place 
strong lower bounds on all of the right-handed squarks. Meanwhile, light EW 
superpartners are generic throughout much of the parameter space. This is espe¬ 
cially the case at lower messenger scales, where a positive threshold correction to 
rrih coming from light Higgsinos and winos is essential in order to satisfy the Higgs 
mass constraint. 
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1 Introduction 

The recent discovery of a Higgs boson near 125 GeV [1, 2] has important and far-reaching 
implications for supersymmetry. In minimal implementations of SUSY (i.e. the MSSM), the 
stops must now either be very heavy (> 10 TeV) or have a large trilinear coupling to the 
Higgs, a so-called ‘A-term' [3-11]. Although the heavy stop scenario is trivial to achieve, it is 
less interesting from both the experimental and the theoretical point of view. Meanwhile the 
large A-terrii scenario allows for stops to be observed at the LHC, and it presents interesting 
challenges for model building. 

While many ideas have been explored on how to generate large, multi-TeV A-terms from 
integrating out the messengers of SUSY-breaking [12-30], perhaps the simplest mechanism 
comes from the MSSM itself - radiatively generating A-terms through the MSSM RGEs. This 
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is especially necessary in the context of gauge mediated supersymmetry breaking (GMSB). 
Indeed, while GMSB naturally solves the SUSY flavor problem and remains one of the most 
well-motivated frameworks for the origin of SUSY breaking at the weak scale (for a review 
and many original references, see [31]), it predicts that the Y-terms are essentially zero at the 
messenger scale. 

In this paper, we will perform a systematic and thorough investigation of GMSB in the pres¬ 
ence of the mh = 125 GeV Higgs constraint. For this purpose, we will employ the framework 
of “General Gauge Mediation” (GGM) developed in [32, 33]. There the model-independent 
parameter space and predictions of gauge mediation were shown to be: 1 

• Flavor universality 

• Negligible Y-terms and B^ 

• (jl is “set by hand” 

• The sfermion soft masses obey the following relations 

2 2 2 
m Hu = m H d = m L 

mq — 2rrifj + m 2 D — m 2 L + ni 2 E = 0 (1.1) 

2 tuq — rn(j — rn 2 D — 2 m 2 L + m 2 E = 0 

All of these conditions hold at the messenger scale M mess and are generally modified by the 
RG-running to the weak scale. They allow for seven independent UV parameters that span 
the full parameter space plus M mess itself, which sets the length of the RG-flow. A convenient 
choice of parameters is 2 

Mi, M 2 , M 3 , 7 Bq, my, m\, and/r (1.2) 

It was shown in [33] that the full GGM parameter space can be realized in terms of weakly 
coupled messenger models. It was further shown in [7] that if one starts with zero Y-terms in 
the UV, then very large gluino masses and high messenger scales are required to generate large 
Y-terms at the weak scale. Our goal in this paper is to build on these works, by exploring 
how the full set of constraints in the UV (the GGM boundary conditions) and the IR (the 
Higgs mass, EWSB and a tachyon-free spectrum) impact the allowed parameter space. In a 
companion paper [35] we will study the corresponding LHC phenomenology. 

1 A common extension of gauge mediation is to include additional Higgs-messenger couplings in order 
to generate /r and B see [34, 19] for a discussion in the context of GGM. This may also generate Y- 
terms and modify the boundary conditions for m 2 H and m 2 Hd . Such models are beyond the scope of 
this work; see section 5.3 for further comments. 

2 In this paper, we assume real gaugino masses and /i to avoid problems with CP violation, but allow 
for both positive and negative values for all the soft masses (including /r) in (1.2). We also assume 
messenger parity in the hidden sector so that U(1)y D-tadpoles are zero in the UV [32], 
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While the GGM parameter space (1.2) is a huge reduction in complexity compared to the 
full 100+ parameters of the MSSM soft SUSY-breaking Lagrangian, it is still challenging to 
survey it fully. The main reason is that the GGM boundary conditions that lead to (1.2) 
are defined in the UV at M mess while all of the other constraints are applied in the IR at 
the weak scale. Even after RG evolving from M mess down to the weak scale and imposing the 
EWSB and Higgs mass conditions, a four dimensional parameter space remains. This is further 
subject to the requirement of a viable (i.e. non-tachyonic) spectrum. Previous attempts have 
dealt with this challenge primarily by taking various 2D slices of the UV parameter space [36— 
40]. Aside from introducing artificial relations among the parameters, this is also suboptimal 
because scanning the GGM parameter space in terms of the UV parameters is in general 
quite inefficient. For instance, the IR constraints might not be automatically satisfied at a 
generic point in the UV parameter space, or the UV parameters might map to uninteresting 
IR parameters, e.g. where some sfermions are extremely heavy and out of reach of the LHC. 

A key idea of this paper is to work directly in terms of an equivalent set of IR soft parameters 
defined at the weak scale: 


Mi, M 2 , A t , rriQ 3 , mfj 3 , m\ 3 , and fj, (1.3) 

In order to efficiently map UV to IR parameters, we make use of a “transfer matrix” approach 
to the MSSM RGEs: for fixed tan /?, M mess and Ms = y JmQ 3 mu 3 , we integrate the RGEs once 
and for all and encode the result as the coefficients of a (bi)linear transformation between UV 
and IR parameters. This approach is quite common in high-scale mediation scenarios, but less 
so in gauge mediation scenarios. Via the transfer matrix and the GGM boundary conditions, all 
other IR parameters are determined in terms of those in (1.3) by a set of algebraic relations. (A 
subset of these relations - those that are one-loop RG invariants - was previously presented and 
studied in detail in [41, 42].) Using these IR relations to reduce the MSSM soft masses to (1.3) 
will streamline the task of scanning over the parameter space, elucidate the phenomenology of 
GGM, and clarify the interplay between the various IR constraints. 

We will take two complementary approaches to exploring the GGM parameter space at the 
weak scale. Our first method is to perform a high-resolution numerical scan on the parameter 
space in (1.3). Since the RGEs depend on M\ only through the small hypercharge coupling, M\ 
plays very little role in the analysis, and so we set M\ = 1 TeV throughout. We explore the role 
of M mess by defining three benchmark scenarios with “low”, “medium” and “high” messenger 
scales, M mess = 10 7 , 10 11 and 10 15 GeV respectively. (Messenger scales higher than 10 15 GeV 
are not considered because gravity-mediated effects are expected to become important, spoiling 
the flavor-universal GGM boundary conditions.) Finally, we choose (mQ 3 ,mfj 3 , M 2 ) to scan 
finely over. For each choice of these parameters, we use the Higgs mass and EWSB conditions 
to eliminate A t , m| 3 and ^ i . We use SoftSUSY [43] to take into account all relevant IR threshold 
corrections. For each point in the stop mass plane, the allowed parameter space is an interval 
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(or collection of intervals) in M 2 . 

We will also study the GGM parameter space analytically in a simplified approximation, 
in order to gain deeper insights. First, we will neglect all of the threshold corrections to 
the EWSB equations and truncate them to tree-level. As we will see, this approximation is 
surprisingly effective. Second, we will greatly simplify the IR relations by using the one-loop 
RGEs and by neglecting contributions from hypercharge and the bottom and tau Yukawas. 
For example, the IR relation for the right-handed slepton mass becomes: 



A complete list of simplified IR relations and their derivation is given in section 2.1. They 
are a central result of this paper, and they will prove to be quite powerful. Together with the 
tree-level EWSB conditions and an accurate Higgs mass calculation via SoftSUSY, we find that 
we can understand nearly all of the features of the GGM parameter space in this approximate 
analytical approach. 

As we will see, the IR relations imply certain orderings of the soft masses. For example, 
we will show that the first and second generation Q and U squarks are always heavier than 
their third generation counterparts, and that the D squarks are always heavier than the lightest 
stop. For the other sparticles, the ordering generally depends on where we are in the parameter 
space. Most importantly, if rriQ 3 < mu 3 , then according to (1.4), left-handed sleptons are 
always lighter than right-handed sleptons. Meanwhile for ?bq 3 > mjj 3 , right-handed sleptons 
are always lighter provided p, is not too large. Based on these orderings, we show that the 
boundaries of GGM parameter space are solely determined by the Higgs mass, EWSB, slepton 
tachyons, and left-handed stop/sbottom tachyons. All other potential constraints (such as 
tachyons from the other scalars) are irrelevant. 

Not surprisingly, the Higgs mass constraint plays an especially important role. The reason 
is that, as noted above, large radiative A-terms with light stops require very heavy gluinos 
in GGM. Such heavy gluinos have a number of effects on other soft parameters through the 
RGEs. For example, as was noted in [7], the stops must be tachyonic at the messenger scale 
and over much of the RG. (See also the nice discussion in the earlier work of [44] and its possible 
implications for fine-tuning.) This is in tension with EWSB, since negative soft masses for the 
stops drive upwards in the RG-running, 

167r2 ^t m H u = 6 Vt( m Q3 + m u 3 ) + ■■■ ( L5 ) 

while EWSB requires < 0 at the weak scale. Of course, the simple way out is to start 
with sufficiently negative m 2 Hu already at the messenger scale. But in models of GMSB, 
m\ 3 = m 2 Hu at the messenger scale, and so left-handed slepton tachyons come into play, 
ruling out combinations of stop masses and A-terms which would otherwise have satisfied the 
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Figure 1: Examples of spectra where both itiq 3 and mjj 3 are small at the weak scale. rriH u is driven 
upward in the RG-running and EWSB is not achieved in the IR (left panel) or the left-handed slepton 
remains tachyonic in the IR (right panel). 

nih = 125 GeV constraint. This logic is further illustrated in fig. 1 for an example point with 
low stop masses. 

Given the role of the stop masses in determining not only the Higgs mass, but also the 
ordering of the slepton masses, the projection of the GGM parameter space into the stop soft- 
mass plane ( m.Q 3 , mjj 3 ) will prove to be extremely useful throughout this paper. A schematic 
representation of the stop mass plane is shown in fig. 2. We have divided it into two halves 
along the diagonal, and we will refer to the mg 3 < rnu 3 ( niQ 3 > mjj 3 ) half as the “LHS” 
(“RHS”) of the stop mass plane. According to our discussion above, on the LHS (RHS), 
left-handed (right-handed) slepton tachyons take precedence in determining the boundaries of 
parameter space. 

These tachyon constraints do not act symmetrically across the diagonal of the stop mass 
plane. In fact, we will show that the right-handed slepton tachyon constraint leads to a 
strict lower bound on mjj 3 of > 1.5 TeV (and becoming even more stringent with decreasing 
messenger scale). Because of the IR relations, there are similarly stringent bounds for all of the 
other right-handed squarks (both up and down-type). Meanwhile, no comparable lower bound 
on m,Q 3 exists on the LHS. Instead, here the boundary arises because a large hierarchy between 
rriQ 3 and M 3 induces a large, negative threshold correction to the left-handed stop/sbottom 
mass, driving it tachyonic. The left-handed squark masses of the first/second generation track 
rriQ 3 . again because of the IR relations. As a result, we find that all three generations of 
left-handed squarks can be arbitrarily light, despite the constraints on GGM parameter space. 

Near the boundaries of parameter space, a convergence of constraints leads to a highly pre¬ 
dictive set of spectra with definite implications for the collider phenomenology. In particular, 
since the boundaries are always determined in part by a sparticle mass going tachyonic, the 
spectrum there is always characterized by light sparticles. For instance, the L and E bound¬ 
aries are always accompanied by relatively light left and right handed sleptons respectively, 
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Figure 2: The anatomy of the IR stop mass plane in GGM. The red line schematically indicates the 
boundary of the viable parameter space; for the actual bounds we refer to fig. 6. The GGM sum rules 
divide the plane along the diagonal, inducing L tachyons on the left and E tachyons on the right of the 
diagonal. On the far left, the parameter space is cut off by left-handed stop/sbottom tachyons coming 
from gluino-induced threshold corrections. 

while the Q boundary predicts light left-handed stops/sbottoms. The tension with EWSB 
generally implies light Higgsinos as well. These are interesting predictions for the boundary 
of the GGM parameter space and provide additional motivation for LHC searches focused on 
stops, sbottoms and EW superpartners. 

Finally, we will see from both the full numerical scan and the approximate analytic approach 
that the sign of ji is an important discrete choice that affects many qualitative features of the 
GGM spectrum. Requiring that there be no pseudoscalar tachyons correlates the sign of /u and 
the range of M 2 throughout GGM parameter space; in particular, only for fi> 0 can one obtain 
M 2 = 0. The Higgs mass receives a ~ 2-3 GeV boost from light charginos and neutralinos 
in the neighborhood of M 2 = 0 , and this leads to a significant decrease in the required At for 
[i > 0 compared to fi < 0. We will see how this difference between the two signs of fj, becomes 
more striking as M mess is decreased and all the various constraints become much stronger. 
Eventually, the positive chargino/neutralino threshold correction becomes critical to fulfill the 
Higgs mass constraint and n > 0 with M 2 ~ 0 dominates the parameter space. Light winos are 
therefore another robust feature of the GGM parameter space with lower M mess , with possibly 
important consequences for the LHC phenomenology. 

The remainder of this paper is organized as follows. In section 2 we discuss the general 
features of the GGM parameter space at the weak scale, and in particular the impact of 
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imposing the Higgs mass constraint. We derive the approximate IR relations that result from 
the GGM boundary conditions in the UV. We use these to identify the relevant constraints 
(EWSB, Higgs, slepton and left-handed stop/sbottom tachyons) and show how they restrict 
the parameter space. Section 3 contains the methodology and results of our numerical scan 
with SoftSUSY. We present results in the stop mass plane, and also along several benchmark 
slices of parameter space, which serve to further illustrate the features of the parameter space 
and the impact of the various constraints. Section 4 is devoted to a detailed semi-analytic study 
of the parameter space, which we use to validate and elucidate features of the full numerical 
scan. We conclude in section 5 with a summary of our results, a brief preview of the upcoming 
paper [35] on GGM collider phenomenology, and a discussion of future directions. Appendix 
A contains validation plots for the transfer matrix and our numerical scan, while appendix B 
discusses in more detail the chargino/neutralino threshold correction to the Higgs mass. 

2 GGM at the weak scale 

2.1 IR relations in GGM 

As described in the introduction, one of the key ideas of this paper is to work directly in terms 
of the IR soft parameters: 


Mi, M 2 , A t , m 2 Q3 , m 2 U3 , m| 3 , /t (2.1) 

We are able to do this efficiently by using a “transfer matrix” approach to the MSSM RGEs. 
For fixed tan/3, M mess , and Ms = mu 3 > we can integrate the MSSM RGEs once and for 

all and relate the UV soft parameters to the IR ones using a set of simple algebraic relations. 
For instance, for M mess = 10 15 GeV, Ms = 3 TeV and tan (3 = 20 we find 3 

m Q3 ~ 0-9 Wq + 2-3 Mf + 0.3 M 2 + ... 

m[/ 3 « 0.8 rhij + 2.0 M 2 + ... (2.2) 

Af ~ —0.99 M 3 — 0.2 M 2 ~\~ ■ ■ ■ 

where ... refers to contributions with smaller coefficients, and the hatted (unhatted) quantities 
are the UV (IR) parameters. To achieve optimal convergence with the remainder of our 
algorithm, we extracted these coefficients using the full 2-loop RGEs of SoftSUSY-3.5.1 [43]. 4 

3 We will fix tan/3 = 20 everywhere in this paper, as this saturates the tree-level contribution to 
the Higgs mass in the MSSM, without being so large that bottom and tau Yukawa effects cannot be 
neglected. As long as tan/3 remains moderately large, we do not expect our conclusions to change much. 

4 The full set of transfer matrix coefficients used in this paper can be accessed in a accompanying 
Mathematica notebook, which is included in the source of this paper on http://arxiv.org/ . 




Using the transfer matrix, we can algebraically reduce all other IR soft parameters to those 
in (2.1). These IR relations are the low-energy versions of the GGM boundary conditions. The 
full set of IR relations is very complicated and we will not reproduce them here. (For a subset 
of these relations that are renormalization group invariants, independent of the messenger 
scale, see the in-depth discussion in [41, 42].) Rather, in this subsection, we will study the IR 
relations in a simplified approximation that consists of using the one-loop RGEs; neglecting 
y 2 , Ur and g\ corrections; and imposing EWSB at large tan/3: 


m l = ~ 2 ("4 + 1/4) + • • • ^ 

Bf, tan/3 = m 2 Hd - m 2 Hu + ... 

These simplified IR relations will form the basis of our understanding of the GGM parameter 
space. Note that they are independent of the Higgs mass constraint; we will come to that 
in the next subsection. All of the sub-leading corrections (and more) are properly taken into 
account in a full numerical scan using SoftSUSY, to be described in section 3. However, as 
we will see through numerous detailed comparisons with this scan, the approximate treatment 
introduced here manages to capture most of the qualitative and even quantitative features of 
the parameter space. 

We begin with the IR relations for the sfermion masses: 


"4,2 f 

«"4 + ^("4 - "4) 

"4,2 C 

s 4 3 + 3 ("4 - "4) 

"4,2 ? 

a "4 

"4,2,3 ? 

- 2^ m Q3 + m %) - 2 mHu 

"4,2,3 ? 

a 2 "4 - \ m Hu + 2^4 


(2.4) 


These relations are satisfied exactly at the messenger scale due to the GGM boundary condi¬ 
tions. In the IR, they are only violated by small effects proportional to y 2 and y 2 . Working 
in the same approximation, we do not concern ourselves with the small splittings amongst 
the three generations of sleptons and right-handed sbottoms (see [45] for a discussion of the 
slepton splitting in GGM). Notice that these relations are independent of the messenger scale 
and the details of the transfer matrix. Thus they are examples of the renormalization group 
invariants discussed in [41, 42]. 

After imposing the large tan /3 EWSB condition rrijj ~ — y 2 , we reduce the other sfermion 
masses to simple combinations of the IR parameters in (2.1). These IR relations have a number 
of interesting consequences, which we list here: 


• The lst/2nd generation Q and U squarks are always heavier than their 3rd generation 
counterparts. We emphasize that this result is not completely trivial once negative rnass- 
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squareds in the UV are allowed (as is the case in GGM), as these could a priori reverse 
the Yukawa effects in the RGEs that usually drive the third generation squarks lighter. 


• The D squarks are always heavier than the root-mean-squared of the stop masses. 

• The right-handed sleptons are strictly heavier than the left-handed sleptons, provided 
that rriQ 3 < mjj 3 ■ For rriQ 3 > mjj 3 , which is lighter depends on ^ 2 . 

Next we turn to the Higgs sector. Here the IR relations are : 5 

2 2 
m H d ~ rn La 

m 2 H u ~ e (5M 2 + dA t ) 2 + am\ 3 - ml (2.5) 

777.(4 ~ -S/x fan p ~ —g 5M-2 g tan /? 

where we have defined 

SM 2 = M 2 + fA t (2.6) 

and 

ml = b{m 2 Q 3 +mjy -cA 2 t (2.7) 

The first relation in (2.5) is a consequence of the GGM boundary conditions at the messenger 
scale and it is only violated by effects. The second and third relations are derived by 
integrating the MSSM RGEs, dropping subdominant contributions proportional to g\. Unlike 
the previous IR relations, these depend on the messenger scale; see table 1 for benchmark 
values of the coefficients a, b, .... In terms of (2.5), the tree-level EWSB equations (2.3) 
become 

e ( 5M 2 + d A t ) 2 + a m\ + g 2 « ml 

3 ( 2 . 8 ) 

— g 8M 2 g tan /3 ~ m ^ 3 + g 2 

From these IR relations, we learn that 

• An important corollary of the formula for 777(4 in (2.5) is that the sign of g and 8M 2 are 
correlated. Concretely, if g < 0 (g > 0) we must have SM 2 > 0 (SM 2 < 0) to avoid 
pseudoscalar tachyons. 

• In fact, pseudoscalar tachyons are always superseded by positivity of m 2 ^ and g 2 , ac¬ 
cording to the second EWSB condition in (2.8). 

• Also from the second line in (2.8) it is clear that g = 0 is not an independent constraint, 
at least in our current approximation, since it always implies 77i 2 3 = 0. 

5 Note that we are assuming A t = B^ = 0 in the UV for simplicity. In practice there are generally 
small, higher-loop contributions to these quantities in GGM. We have checked that none of our results 
depend sensitively on A t and B M being literally zero in the UV. 
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■Mmess 

a 

b 

c 

d 

e 

f 

g 

P 

q 

10 15 GeV 

0.59 

0.43 

0.46 

0.7 

0.2 

0.73 

0.45 

1.62 

0.47 

10 11 GeV 

0.69 

0.33 

0.49 

0.83 

0.1 

0.87 

0.28 

1.95 

0.39 

10 7 GeV 

0.83 

0.18 

0.44 

1.51 

0.02 

1.09 

0.13 

3.08 

0.31 


Table 1: Parameters used in the IR relations (2.5)-(2.10) for various values of M mess , with tan/3 = 20 
and Ms = 3 TeY. 


• From (2.4) and (2.8), it follows that 


mi 3 < 




(2.9) 


Since | for all messenger scales, we expect that E tachyons are always a stronger 
constraint than C /3 tachyons. Ultimately this translates into a strong lower bound on 
rriu 3 , as we will show in section 3. 


• The quantity ttiq defined in (2.7) must be positive, otherwise the first EWSB condition 
in (2.8) cannot be satisfied with non-tachyonic sleptons. This places an upper bound on 
the magnitude of the A-term allowed at each point in the stop mass plane. 


Finally, let us comment on the role played by the gluino. The IR gluino mass is given in 
terms of At and SM 2 by 

M 3 « -(pA t + qSM 2 ) ( 2 -10) 

where benchmark values of p and q are listed in table 1 . This equation shows how M 3 is linearly 
related to At and 5M 2 . The proportionality constant p moreover increases with lowered M mess . 
This reflects the fact that a larger gluino mass is needed to achieve the same At for a shorter 
amount of RG running. As we will see in the following sections, enormous gluino masses are 
generally required to achieve the large A-term scenario with lower messenger scales, and this 
can result in large gluino-induced threshold corrections to the IR squark masses, as given by 
equation (34) in [46]: 

<5m| = ^p^m| ^1 + 3x + (x — l ) 2 log \x — 1| — x 2 logic + 2xlog 

Here rriq stands for any of the squark soft masses, and x = M|/m|. These threshold corrections 
are generally negative for the gluino masses of interest (i.e for M 3 Mg), and will eventually 
turn the physical squark mass tachyonic. As we will see, this effect is ultimately responsible 
for the left-most boundary in fig. 2 . 


Mf 

9 

m~ 


( 2 . 11 ) 
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2.2 Imposing the Higgs mass constraint 


Now we will impose the Higgs mass constraint and discuss its implications for GGM. Through¬ 
out this work, we will require = 123 GeV, in order to account conservatively for the theory 
uncertainty [47] in the Higgs mass calculation. In the MSSM, the Higgs mass is given by the 
well-known formula 

3v 2 ( 

m h = m| cos 2 (2/3) + ^2 ( M 4 l°g 


Ms 


m+ 


A? 


-4 +T7, \yt\~ 


Mjg 



+ ... 


( 2 . 12 ) 


Here the ... denote important additional corrections from mg 3 ^ mu 3 -, other sparticle thresh¬ 
olds and higher loops. These are accounted for in our analysis using SoftSUSY. 

The Higgs mass stringently constrains the stop masses and the A-terms in the MSSM; for 
TeV-scale stops, the A-terms must generally be multi-TeV. An example of this is given in fig. 3. 
Shown here are contours of the “normalized A-term” 


Rt = 


\M 

yJ m Q 3 + m u 3 


(2.13) 


required for mh = 123 GeV in SoftSUSY, with all other superpartner masses set to Mg. (l As 
we lower the stop masses, the required Rt increases, and for stops below ~ 1 TeV, the Higgs 
mass constraint cannot be satisfied. 

In GGM the Higgs mass constraint is even more stringent, since we cannot vary the other 
parameters of the MSSM arbitrarily. As noted in the previous subsection, for m q < 0, the first 
EWSB equation in (2.8) has no valid solution, since the LHS must be a sum of non-negative 
quantities. This translates to the requirement that 


R% < b/c (2.14) 

and from table 1, we see that \Jbf~c = 1.01, 0.85, 0.69 for M mess = 10 15 , 10 11 , 10 6 7 * GeV respec¬ 
tively. Comparing with fig. 3, we see that stop masses that would otherwise be allowed by the 
Higgs mass constraint are ruled out in GGM by the combination of EWSB and no-tachyon 
conditions. (Keep in mind that while m ( 2 = 0 furnishes an absolute boundary to the parameter 
space, there can be even more stringent boundaries due to tachyon constraints.) 

Clearly, the decrease of \Jb/c with messenger scale amplifies the tension between large 
A-terms and the other constraints. This will serve to enhance the role of secondary threshold 
corrections that can increase mh and allow for smaller At. As we will see in section 3.3 (and will 
discuss further in appendix B), chief among these is a ~ 2 GeV positive threshold correction 

6 We emphasize that this figure is meant to give a general impression and should not be taken literally. 

The A-term required for mh = 123 GeV can depend sensitively on the masses of the other superpartners 

and their contributions to the Higgs mass. 
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m Q3 (TeV) 


Figure 3: Contour plot of the R t needed in the MSSM for rrih = 123 GeV as a function of the stop 
soft masses. Other soft masses were set to Ms = 3 JuiQ 3 mu 3 - The gray regions indicate places where 
the Higgs mass constraint was impossible to satisfy. In the lower left corner, this is because no A-term 
satisfies the constraint. On the sides, this is because various threshold corrections coming from the very 
split stops render the output of SoftSUSY unreliable. 

to uih coming from light winos and Higgsinos when M 2 and /r are both close to zero. Since 
the M 2 ~ 0 region requires (jl > 0 according to the first bullet point below (2.8), this will lead 
to a marked difference between fi < 0 and fi > 0 parameter spaces as M mess decreases. 

To summarize, we have seen in this section that the IR soft parameters of GGM are 
related to those in (2.1) via a set of simple algebraic relations. Some of these IR relations are 
renormalization group invariants along the lines of [41, 42], while others are not. Using these 
relations, we have shown that the IR soft masses obey certain fixed orderings. In particular, 
the only soft masses in GGM that can become tachyonic independently of others are itiq 3 , 
and m| 3 . Ah other soft masses are always positive as long as these soft masses are positive. 
Together with the Higgs mass and EWSB constraints, these determine the boundaries of GGM 
parameter space. 

We also showed that important qualitative dividing lines cutting through the parameter 
space include: the diagonal of the stop mass plane, the sign of / 1 , and the (anti-correlated) 
sign of SM 2 . Using the important variables m q and Rt, we saw how decreasing the messenger 
scale results in increasingly stringent constraints on the parameter space. In the following 
sections, we will confirm this general picture using a high-resolution numerical scan of the 
GGM parameter space, together with an analytical approach based on the approximate IR 
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relations and tree-level EWSB conditions. 


3 Scanning the GGM parameter space 

3.1 Details of the scan 

In the introduction, we sketched out the steps in our numerical scan of the GGM parameter 
space. These steps are summarized in fig. 4. Here, we will describe them in more detail. 

1. We define “low”, “medium” and “high” messenger scale benchmarks corresponding to 
M mess = 10', 10 11 and 10 15 GeV respectively. Moreover, having verified that a\ effects 
have very little impact on the analysis we set Mi = 1 TeV throughout this paper. 

2. We choose to eliminate A t , m\ 3 and /x using the Higgs mass and EWSB equations, since 
the former depends strongly on At, while the latter are sensitive to and //. An 
additional benefit of this choice is that appears linearly in the EWSB equations 
( 2 . 8 ). 

3. This leaves niQ 3 , mjj 3 and M 2 as independent parameters. As described in the introduc¬ 
tion, a convenient way to view this remaining parameter space is that for every point in 
the stop mass plane, all soft parameters are functions defined on an interval or collection 
of intervals in M 2 • The A 0 tachyon condition cuts the M 2 interval into two disconnected 
pieces, one for each sign of /x. Both pieces are further bounded by requiring the absence 
of slepton tachyons and by the fi 2 > 0 conditions. These features are illustrated in fig. 5 
for an example point in the stop mass plane. 

4. Finally, we perform a high-resolution three-dimensional scan over (mq 3 , rriu 3 , M 2 ). Near 
the boundary of the parameter space the resolution of the scan is further increased, such 




Mi m Qz ,m U3 ,M 2 


m L3 , M 


10 7 ,10 n ,10 15 

GeV 


1 TeV fine scan mu = 123 GeV EWSB conditions 


Figure 4: A schematic diagram illustrating the steps in our numerical scan over the GGM parameter 
space. We trade all UV parameters for the more physical IR parameters in (2.1) using the transfer 
matrix and use the EWSB conditions and the Higgs mass constraint to eliminate A t , m? L and /x. We 
further fix M\ = 1 TeV, tan/3 = 20, scan coarsely over M mess , and scan finely over toq 3 , mjj 3 and M 2 . 
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Figure 5: Example of the behavior of various soft parameters over the M 2 interval. Here toq 3 = 1 
TeV, mu 3 = 3 TeV and M mess = 10 15 GeV. The feature around M 2 ~ 0 is due to threshold correction 
to rrih and will be discussed in detail in the next sections. 

that this important region is sampled as accurately as possible. The end result is a 
complete grid of valid spectra spanning the GGM parameter space and satisfying the 
Higgs mass and all other IR constraints. Appendix A contains several validation plots 
which demonstrate the convergence of our algorithm. 

In practice, step 2 above is the most challenging part of the analysis, because of the 
complicated threshold corrections that must be taken into account. Because we use SoftSUSY 
to implement all the threshold corrections and RGEs, there is a particular order in which 
we must solve these constraints. For a given input and at the messenger scale, 

SoftSUSY imposes the EWSB conditions and returns fj, and B M at the messenger scale. It also 
computes mh • Thus for each (mQ 3 ,mjj 3 , M 2 ), fi is given automatically by SoftSUSY, but we 
must perform an auxiliary 2D scan over At and and numerically solve the -B^(UV) = 0 
and rrih = 123 GeV constraints. In principle, this adds two extra dimensions to our scan. A 
brute force, flat scan over (A t , rrij 3 ) proves to be computationally unfeasible. Instead, we use 
an iterative method where we sample a few points in the neighborhood of a seed guess, linearly 
interpolate in B fl ( UV) and mh to establish a new seed, and repeat. Typically this converges 
quickly, after just a few steps, on an extremely accurate solution. Of course, for this to work, it 
is crucial to obtain an accurate initial seed for (At, ). We accomplish this by extrapolating 
from a neighboring point in the (mQ 3 ,mjj 3 , M 2 ) parameter space. 
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3.2 Results: a “birds-eye view” 


In the rest of this section, we will exhibit the results of the scan outlined above. The primary 
focus here will be on describing its features; a fuller analytic understanding in terms of the IR 
relations (2.4) and (2.5) and the tree-level EWSB equations (2.8) will follow in section 4. 

We begin with a “birds-eye view” of the parameter space: the viable region projected onto 
the stop mass plane. Shown in fig. 6 is the full result of our numerical scan for the three 
different values of M mess and the two signs of /i. There are several interesting features of these 
plots which highlight the general points made in section 2. These include: 

• As expected from the discussion in section 2.2, the allowed region shrinks as M mess 
decreases. 

• For M mess = 10 15 GeV, the difference between fi < 0 (blue) and /r > 0 (orange) is 
minimal, but it becomes increasingly dramatic as M mess decreases. As we will see in 
more detail below, this is due to the increasing importance of the chargino/neutralino 
threshold correction to m^. 

• Although the lower bound on rag 3 becomes increasingly stringent with lower M mess , the 
physical mass of the mostly-left-handed stop can be arbitrarily low. We will verify in 
section 3.4 that this is due to the gluino threshold correction. 

• The same is not true for the mostly-right-handed stop, whose mass is bounded from 
below by ~ (1.5, 2, 2.5) TeV for M mess = (10 15 , 10 11 , 10 7 ) GeV. According to (2.9), the 
right-handed slepton tachyon constraint prevents mu 3 from becoming too light. 

In the remainder of this section, we will further elaborate on these and other features by 
“zooming in” on these plots and exploring the parameter space along three different benchmark 
lines depicted in fig. 7. These lines are chosen in order to illustrate the behavior of the 
parameter space as we approach the L, E and Q tachyon boundaries. Since the allowed 
parameter space for M mess = 10 7 GeV is smaller, we will focus on the M mess = 10 15 GeV and 
Mmess = 10 11 GeV cases. 

3.3 The L and E lines 

These lines are defined by varying mu 3 while holding mg 3 fixed to the benchmark values shown 
in fig. 7. The physics along the L and E lines is qualitatively similar, so we will consider both 
simultaneously in this subsection. 

Shown in fig. 8 are plots of the normalized A-term Rt vs. mu 3 for these lines in the stop mass 
plane. As mjj 3 approaches the boundary (i.e. its minimal allowed value), there are two features 
worth noticing: first, Rt increases due to the Higgs mass constraint, as expected from fig. 3. 
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m Q3 (TeV) 

Figure 6: Allowed regions in the stop soft mass plane for /i < 0 (/r > 0) indicated by the blue (orange) 
shaded regions. The gray dots are the allowed physical stop masses, which can differ significantly from 
the soft masses due to the gluino threshold correction. The wedge along the diagonal is a result of the 
level repulsion between the two stop mass eigenstates. 


This leads to an increasing tension with EWSB, as explained in section 2.2. Second, there is a 
range of Rt values for each mjj 3 , which occurs because we marginalized over the M 2 interval in 
this figure. The range for Rt shrinks to zero once mu 3 approaches its minimal allowed value. 
This indicates that the viable M 2 interval shrinks to a point prior to disappearing completely. 
Figs. 8 and 9 also illustrate very starkly the difference between jx < 0 and jx > 0: we see 
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m Q 3 = 1 TeV 
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10 11 GeV 

771q 3 = 1.5 TeV 
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m Q3 = 4 TeV 
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10 11 GeV 

m Q3 = 4 TeV 

Q 

10 15 GeV 

mjj 3 = 4 TeV 

Q 

10 11 GeV 

mjj 3 = 4 TeV 


Figure 7: The same schematic cartoon of the GGM parameter space shown in fig. 2, now with three 
lines labeled by ( Q ), (L) and (E) indicating the three benchmark lines we study in detail in this section 
and the next. 


that the A-terms are mostly constant across the M 2 interval for /j, < 0, as would be expected 
from the Higgs mass constraint, but for /i > 0 they vary quite a lot across the M 2 interval. 
Evidently, the magnitude of the A-term required for nif l = 123 GeV decreases significantly in 
the neighborhood of M 2 = 0. 

The reason for this decrease can be traced back to a positive one-loop threshold correction 
to the Higgs mass coming from loops of light Higgsinos and winos, see appendix B for more 
details. As the magnitude of M 2 decreases, rrih grows by ~ 2-3 GeV, and this greatly relaxes 
the demands on the A-term. Since M 2 = 0 is only accessible for fj, > 0 due to the pseudoscalar 
tachyon constraint (first bullet point below (2.8)), only fi > 0 is sensitive to this threshold 
correction. This explains why the allowed parameter space (fig. 6) for ^ > 0 becomes much 
larger than the one for fj, < 0 as M mess decreases. As the constraints on the stop mass plane 
become more stringent, the importance of the small-M 2 threshold correction is magnified. To 
the point that for M mess = 10' GeV, the constraints basically kill off the entire parameter 
space, except where this small-M 2 threshold correction is present. 

Finally, in figs. 10 and 11 are plots of p, and the relevant slepton mass across the M 2 
interval, again with mjj 3 varying along the L or E line as indicated by the color coding. In 
these figures the correlation between the sign of /r and the viable range of M 2 is especially 
manifest. We highlight some other general features of these plots. For /r < 0: 

• \n\ is always monotonically decreasing with M 2 . 

• The lower end of the M 2 interval is determined by —> 0 or m,E 3 —> 0 on the L or E 
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line respectively, always with /z / 0. 

• The upper end is determined by /j. —>• 0, and on the E line it is sometimes accompanied 
by mE 3 —> 0 (i.e. for M mess = 10 15 GeV and lowering mu 3 closer to the boundary). 

Meanwhile, for ju > 0: 

• Again, the M 2 = 0 region has a large effect on the plots. |/z| is no longer monotonic but 
tends to rise and fall as we cross M 2 = 0. 

• In all cases, the allowed M 2 interval starts to center around M 2 = 0 as mu 3 is lowered. 
When this happens, the M 2 interval becomes bounded by // —> 0 and rriE 3 —> 0 on both 
ends along L and E lines respectively. 

In section 4, we will understand these features analytically in terms of the approximate IR 
relations and tree-level EWSB equations described in section 2.1. 
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Rt M mess = 10 15 GeV 



m U3 (TeV) 


Rt M mess = 10 11 GeV 
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Figure 8: R t as a function of mjj 3 along the L (top) and E (bottom) benchmark lines, for M mess = 
10 15 GeV (left) and 10 11 GeV (right). The orange and blue shaded points correspond to the /r > 0 and 
/i < 0 branches respectively. 
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1.1 




Figure 9: Plots of Rt vs. M 2 along the L (top) and E (bottom) benchmark lines, for M mesa = 10 15 GeV 
(left) and 10 11 GeV (right). The orange and blue shaded points correspond to the fi > 0 and /x < 0 
branches respectively. The progressively darker shading of the colors indicates decreasing mu 3 ■ 
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Figure 10: Plots of |/i| (top), itil 3 (bottom) vs. M 2 along the L benchmark line, for M mess = 10 15 GeV 
(left) and 10 11 GeV (right). Color schemes are as in fig. 9. The black curves correspond to the 
quantitative predictions of the lower end of the M 2 interval for /i < 0 in the semi-analytic approximation, 
see section 4 for details. 
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Figure 11: Plots of \/i\ (top), itie 3 (bottom) vs. M 2 along the E benchmark line, for M mess = 10 15 GeV 
(left) and 10 11 GeV (right). Color schemes are as in fig. 9. The black curves correspond to the 
quantitative predictions of the lower end of the M 2 interval for p < 0 in the semi-analytic approximation, 
see section 4 for details. 
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3.4 The Q line 


Finally, we come to the Q benchmark line shown in fig. 7. Here we fix mu 3 = 4 TeV for both 
M-mess = 10 15 GeV and M mess = 10 11 GeV. The normalized A-terrn, /r and slepton masses 
all exhibit the same general behavior along the M 2 interval as on the L line, so we will not 
show these plots again for the Q line. The big difference with the L and E lines is that the 
Q line is not cut off by EWSB and slepton tachyons, but rather by a Q tachyon. Shown in 
fig. 12 is the pole mass of the lightest stop vs rriQ 3 along the Q line. We see that begins 
to differ significantly from the soft mass niQ 3 as we approach the boundary of the stop mass 
plane, ultimately decreasing to zero. (A similar effect occurs for the pole mass of m^.) As in 
fig. 6 , we see that the mostly left-handed stop mass eigenstate can be arbitrarily light despite 
the Higgs mass constraint. 

Also shown in fig. 12 is the range across the M 2 interval of mg 3 subject to the gluino 
threshold correction (2.11). We see that it agrees quite well with the full numerical result 
given by SoftSUSY. This confirms that the gluino loops dominate the threshold corrections to 
the lightest stop mass and are ultimately responsible for rri^ going tachyonic at low rag 3 . 

The plots in fig. 13 illustrate more properties of the gluino mass along the Q line. For 
fj, < 0 where the A-term is mostly constant, we see the tight linear relationship between M 2 
and M 3 encoded in equation (2.10). Furthermore, we see that the Higgs mass constraint forces 


m-{ TeV) M M ess=10 15 GeV 



m-(TeV) M Me ss=10 11 GeV 



Figure 12: Plots of the physical stop pole mass as a function of itiq 3 along the Q benchmark for 
/i > 0. (The [i < 0 plots show similar behavior.) The blue shaded region shows the variation within 
the allowed M 2 interval of itiq 3 subject to the gluino threshold correction (2.11). 
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Figure 13: Plots of M 3 vs. M 2 along the Q benchmark line, for M mess = 10 15 GeV (left) and 10 11 GeV 
(right). Color schemes are as in fig. 9. 


the gluino mass to be quite large overall, and causes it to grow as rriQ 3 is lowered. Also the 
gluino mass is generally larger for M mess = 10 11 GeV than for M mess = 10 15 GeV, since a 
larger M 3 is needed to obtain the desired A-term as predicted by equation (2.10). The large 
hierarchy between rag 3 and M 3 enhances the gluino threshold correction for low values of m. q 3 
and lower messenger scales. This is the reason for the increasing lower bound on mg 3 in the 
plots in fig. 6 . 


4 Interpretation 

In this section, we will understand the features of the numerical solution in terms of the tree- 
level EWSB equations (2.8) and the IR relations (2.4) and (2.5) . Having achieved an analytical 
understanding of the Q tachyon boundary through ( 2 . 11 ) in the previous subsection, we will 
focus on how the EWSB and slepton tachyon boundaries are determined analytically. 

We will organize our discussion in this section around the sign of /r. In previous sections, 
we have seen repeatedly that the (i < 0 and > 0 branches differ qualitatively due to the 
presence of the M 2 = 0 threshold correction to the Higgs mass for the latter branch. As a 
consequence, At is basically constant across the M 2 interval for n < 0 , while this is not the 
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case for the g > 0 branch. For g < 0, this gives us much greater control over the parameter 
space, since we can fix At in all the equations above. The g > 0 branch on the other hand 
requires greater care, and we will consider it separately. 

For the convenience of the reader, we repeat here the most important formulas and def¬ 
initions from section 2.2 used in the analysis below. The EWSB equations can be written 
as 

e (SM -2 + d A t ) 2 + a rnj + g 2 « 

3 (4.1) 

— g 5M 2 g tan /3 « + g 2 

with 


5M 2 = M 2 + fA t 
ml = b (■ m 2 Q 3 + mfj 3 ) -cA 2 t . 

The {a, 6, c, d , e, /, g} are the numerical constants in table 1 and are determined by the RGEs. 
The soft mass of the right-handed slepton is furthermore given by 

™l 3 ~2m| 3 + ^/x 2 + ^(m^3 -m 2 Q3 ). (4.3) 

In the remainder of this section, we will heavily rely on these relations. We further make 
the following approximations: all of our formulas in this section will be to first non-trivial 
order in the 1/tan/3 expansion; we are using the tree-level EWSB equations; we are neglecting 
corrections proportional to gi, yi, and y T . Also, for the most part, we will ignore the mild 
variation of the {a, b, c, d, e, /, g} coefficients across the stop mass plane. 


4.1 fi < 0: characterizing the M 2 interval 


The analysis presented in this subsection and the next applies to points with g < 0, which are 
the blue shaded points in the plots in section 3.2. We begin our discussion in the bulk of the 
stop mass plane, i.e. with large mg. Here it is convenient to introduce a new variable: 

m 2 = m§ - ^a(m 2 Q3 - m 2 U3 )9(m 2 Q3 - mfj 3 ) (4.4) 

where 9 is the Heaviside step function and m-Q is defined in (4.2). So m 2 = m q on the LHS of 
the stop mass plane, while it equals ?u,q — |a(r?ig 3 — mjj 3 ) on the RHS of the stop mass plane. 

On the LHS of the stop mass plane, we expect the M 2 interval is bounded by left-handed 
slepton tachyons. Setting m^ 3 = 0 on the LHS, we find two solutions to (4.1), one at small 
6M 2 : 


5M 2 = —-—- \/m 2 — ed 2 A? 
g tan p v 

g = —Jm 2 — ed 2 A 2 


(4.5) 
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and one at large 5M 2 '. 


5M 2 


—d At + 



g = 0 


(4.6) 


The solutions in (4.5) and (4.6) should thus correspond to the two endpoints of the M 2 interval. 

Sufficiently far into the RHS of the stop mass plane, the M 2 interval should be bounded by 
right-handed slepton tachyons. Setting m 2 E3 = 0 on the RHS again yields two solutions, one 
at small <5M 2 : 


SM 2 


A* = - 


3 m? — e d 2 A 2 + a'(mg 3 — m^ 3 ) 

4 g tan fiyj a! ( m 2 — e d 2 A |) 

j m 2 — e d 2 A 2 

' a’ 


and one at large 5M 2 : 


5M 2 = -dA t + 

= 3 m? Q3 - mfj 3 

4 g 5M 2 tan (3 



(4.7) 


(4.8) 


with a! = 1 — a/4. 

In general, the approximate solutions (4.5) and (4.7) correctly characterize the lower end¬ 
point of the M 2 interval and the general trends along the M 2 interval, but (4.6) and (4.8) fail 
to characterize the behavior at the upper endpoint of the M 2 interval. In more detail: 


• According to the approximate solutions, the lower endpoint of the M 2 interval is char¬ 
acterized by m | 3 —> 0 or m 2 E3 —> 0 with /i 7 ^ 0. These features are all clearly borne 
out in the full solution, see figs. 10 and 11. In these figures, we have also indicated the 
quantitative predictions of (4.5) and (4.7) for the location of the lower M 2 endpoint, as 
shown by the black line. We see that it describes the full solution well. 

• As 5 M 2 increases, it is straightforward to show from the EWSB equations (4.1) that \/i\ 
always monotonically decreases with 5M 2 , while m / 3 and m 2 E3 must rise and fall. These 
trends are clearly borne out in figs. 10 and 11 . 

• Finally, at the upper endpoint of the M 2 interval, the approximate solutions in (4.6) and 
(4.8) predict m 2 L —> 0 or m 2 E3 —> 0 with tan ^-suppressed (or zero) ji. While g does 
become quite small in general, we observe that m\ 3 and m 2 Es appear to be cut off at a 
large value at the upper end of the M 2 interval. This can be traced back to the first 
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EWSB relation in (4.1) which implies 


m ^ 3 ~ —g 5M-2 g tan j3 

3 

m 2 E3 ~ -2g SM 2 g tan (3 - -(m,Q 3 - mjj 3 ) 


(4.9) 


when g is small. The factor of tan (5 and large 5M 2 mean that g has to become extremely 
small before mx , 3 and rn e 3 start to visibly approach zero. Such small values of g become 
sensitive to various effects we have neglected, e.g. threshold corrections to the tree- 
level EWSB equations, and the finite-resolution effects of our grid. Evidently, these are 
enough to cut out the ?bl 3 —> 0 and mE 3 —> 0 behavior at large M 2 . 


• Correspondingly, we find that the quantitative predictions for the upper endpoint given 
in (4.6) and (4.8) do not work so well since these assumed = 0 and vr? E = 0 from 
the start. 


4.2 fi < 0: Approaching the boundaries 

Having discussed the behavior of the M 2 interval in the bulk of the stop mass plane, now we 
turn to its behavior as we approach the boundaries of the stop mass plane, i.e. as we decrease 
m 2 . Shown in fig. 14 are contours of m 2 — ed 2 A 2 for M mess = 10 15 and 10 11 GeV.' At the 
zero contour, the approximate solutions (4.5) and (4.7) become imaginary and are no longer 
valid. We see that the zero contour does a fairly good job of characterizing the boundary of 
the stop mass plane. We have verified that the largest discrepancies for M mess = 10 15 GeV 
(LHS) arise due to sub-leading effects that we have neglected in this simplified semi-analytic 
treatment, specifically corrections proportional to Mi and the variation of the transfer matrix 
along the stop mass plane. 

We must address one technicality, however, before declaring victory. For m 2 < ed 2 A 2 , 

the approximate solutions actually have a second phase where the M 2 interval is bounded by 
/2 

5M 2 = —dAt ± i/y. This phase is distinguished by small g throughout the M 2 interval; in 
fact, on the LHS, g goes to zero at both ends and is non-monotonic on the interval. Because of 
the very small g, we expect this entire phase to not be robust against threshold corrections and 
finite-resolution effects. Indeed, we find that the first phase seems to dominate the parameter 
space of the full numerical solution, and we only see any evidence for the second phase in a 
tiny sliver of the LHS of the stop mass plane for M mess = 10 15 GeV. In any event, the question 
as to whether this phase exists or not is mainly academic, since it would be largely excluded 
by the LEP bound on charginos. 

'Again, we don’t show M mess = 10' GeV here because it is a very tiny region for g < 0 and 
appears to be subject to threshold corrections from enormous gluino masses that make the semi-analytic 
approximation unreliable. 
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Figure 14: Viable points in the stop mass plane for g < 0, with contours of constant m 2 — ed 2 A 2 in 
TeV 2 (thin) overlaid. The blue dots are the individual points of our full numerical scan, to illustrate our 
resolution. We see that the boundary of the stop mass plane is characterized fairly well by m 2 = ed 2 A 2 
(thick). 


4.3 n > 0: the role of the M 2 ~ 0 region 


Next we turn to the g > 0 case. With just a constant A-term, the situation for g > 0 would 
be nearly identical to that of g < 0. However, we have seen in figs. 8 and 9 that the A-term 
needed to achieve rrih = 123 GeV depends strongly on M 2 in the neighborhood of M 2 = 0. As 
explained in appendix B, this is due to the influence of one-loop chargino-neutralino threshold 
corrections. Since these threshold corrections are positive, the required A-terrn is significantly 
smaller around M 2 ~ 0. A smaller A-terrn relieves the tension with tachyons and EWSB, and 
so this can allow points to survive further into the stop mass plane than would otherwise be 
the case. We have already seen this illustrated in fig. 6 , most dramatically for M mess = 7. 

Here we will understand this in more detail using the approximate tree-level equations. To 
begin, let’s set M 2 = 0 (i.e. 5M 2 = / A t ) in (4.1). We see that m | 3 is always large and positive 
in this regime. Solving for g and m | 3 we find 


ml - e(d + f) 2 A? 

H{M 2 = 0) = 0 M n + • • • 


agf {-A t ) tan /3 


m\ (M 2 = 0) = 2 


ml - e(d + f) 2 A: 2 - f a(m^ 3 - m^) 


(4.10) 


+ ... 


where ... denote higher order corrections in 1/tan f3. These are monotonically decreasing 
functions of —At. As we move away from M 2 = 0, the value of —At needed to realize m/,, = 
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123 GeV increases significantly, so this has the effect of decreasing /x and . This results in 
the rise and fall of these parameters around M 2 observed in figs. 10 and 11 . 

As we decrease TOq, the values of /x(M 2 = 0) and rn 2 E ^ (M 2 = 0) decrease, until eventually 
they cross zero. Since the A-terrn here is so much smaller than away from M 2 = 0, this can 
occur further into the stop mass plane than the boundaries discussed in the previous subsection. 
In fig. 15 we show the allowed points for /x > 0, with the zero contours of /x(M 2 = 0) and 
r n E 3 (M 2 = 0) overlaid as predicted by the semi-analytic method. (To improve the accuracy of 



Figure 15: Viable points in the stop mass plane for fi > 0, with contours of /x(M 2 = 0) = 0 (dashed) 
and m 2 E;s (M 2 = 0) = 0 (solid) overlaid. The orange dots are the individual points of our full numerical 
scan, to illustrate our resolution. 
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these contours, we have included the sub-leading 1/tan/3 corrections in (4.10).) We see that 
this does an excellent job of characterizing the boundary of the stop mass plane for /r > 0. 

5 Discussion 

5.1 Summary 

In this paper, we have initiated a comprehensive study of the GGM parameter space following 
the discovery of the Higgs at mh = 125 GeV. In pure GGM, we have at the messenger scale: fla¬ 
vor universality, three independent gaugino masses, = At = 0, and fi “set by hand”. Using 
a transfer matrix approach to the MSSM RGEs, we traded the GGM parameter space defined 
at the messenger scale M mess for a set of IR variables (niQ 3 , mjj 3 , mt 3 , Mi, M 2 , A t , /j). The 
RGEs depend only weakly on Mi through the hypercharge coupling, and the results do not 
qualitatively depend on it. We therefore fixed Mi = 1 TeV throughout the analysis. We fur¬ 
thermore chose three benchmark values of M mess = 10 1 ' 5 , 10 11 , 10 7 GeV. Then the EWSB and 
Higgs mass constraints eliminated three variables, leaving us to scan over (mQ 3 , mjj 3 , M 2 ). 
At each point in the stop mass plane, the viable parameter space was an interval in M 2 . 

We performed a detailed numerical scan and validated it analytically using a set of approx¬ 
imate IR relations together with the tree-level EWSB conditions. Through these methods, 
we thoroughly explored the physically viable parameter space, with a special focus on how 
its boundaries are determined. By using the IR relations to express all other IR soft masses 
in terms of (1.3), we showed that the only relevant constraints on the GGM parameter space 
(besides EWSB and the Higgs mass) were slepton tachyons and left-handed stop/sbottom 
tachyons. All other scalar tachyons were less constraining. Furthermore, we showed that for 
tyiq 3 < mu 3 , only left-handed slepton tachyons and left-handed stop/sbottom tachyons were 
relevant, while for ?t?.q 3 > mjj 3 only right-handed slepton tachyons were relevant. Finally, we 
saw that the constraints became more stringent and the viable parameter space smaller as 
Mmess w as lowered. 

One of the most striking results of our analysis was an absolute lower bound on right-handed 
stop mass, coming from a combination of EWSB and slepton tachyon constraints. In particular 
mu 3 > 1.5, 2, 2.5 TeV for M mess = 10 15 , 10 11 , 10 7 GeV respectively. However, for any mes¬ 
senger scale, the left-handed stop and sbottom could be arbitrarily light. The constraints were 
always such that the optimal point for mh = 125 GeV with ?t?.q 3 ~ mu 3 ~ \At\/V& ~ 1 TeV 
could not be achieved in GGM. 

We also identified the sign of (j, as playing an important role in the qualitative behavior of 
the parameter space. For ^ < 0, the A-term was mostly constant across the M 2 interval, being 
determined by the Higgs mass constraint. However, for fj, > 0, a positive threshold correction 
to mh coming from light charginos and neutralinos allowed the A-terrn to be much smaller in 
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the neighborhood of M 2 ~ 0. This played an especially important role for M mess = 10' GeV, 
where the tension between the Higgs mass constraint and EWSB and tachyons was so strong 
that essentially the only viable parameter space had /r > 0 and M 2 ~ 0. 

5.2 Preview of the LHC phenomenology 

In a companion paper [35], we will explore the LHC phenomenology of the GGM parameter 
space with = 125 GeV. In this subsection we will give a brief preview. 

Our semi-analytic understanding of the GGM parameter space gives us great control over 
the spectrum as we move around in the stop mass plane and the M 2 interval. In particular, it 
allows us to understand under which conditions a given SUSY particle can be light. On the one 
hand, this gives us sharp predictions for the properties of the NLSP, and on the other hand 
it singles out the dominant production channels. Together these two pieces of information 
determine most of the collider phenomenology. 

We saw in our analysis how the constraints became increasingly more stringent as Al mess 
was decreased. Assuming the conventional relation between NLSP lifetime and the messenger 
scale (see e.g. [31] for a review), this implies a strong preference for long-lived NLSPs at the 
LHC. While long-lived neutral NLSPs escape the detector without leaving any track, long- 
lived NLSPs carrying SM charges are very well constrained at the LHC by inclusive CHAMP 
searches. These bounds will play a substantial role in constraining the GGM parameter space. 
For the lowest value of M mess that we considered ( M mess = 10' GeV), the NLSP decay to 
the gravitino may be non-prompt but still inside the detector volume. Constraining these 
scenarios is an interesting challenge for LHC searches (see [48] for a recent discussion) and our 
work further motivates efforts to improve coverage at Run II. 

The dominant component of the colored production cross section will come from left- 
handed squarks throughout much of the GGM parameter space. We showed already that the 
left-handed stop/sbottom can be arbitrarily light. The IR relations (2.4) indicate that the 
left-handed squarks of the first and second generations are heavier, but there are points on the 
M 2 interval where /i ~ mi 3 ~ 0, so they can also become very light. This can be verified in our 
full numerical scan. Meanwhile, the gluinos are generally forced to be very heavy by the Higgs 
mass constraint, especially at lower stop masses and/or lower messenger scales. (Of course, 
with sufficiently heavy stops, A-terms are not required for the Higgs mass, and then the gluino 
can be arbitrarily light.) Finally we saw how the right-handed stops are always forced to be 
at least 1.5 TeV due to the right-handed slepton tachyon constraint. The IR relations (2.4) 
imply that the lst/2nd generation right-handed up squarks are even heavier. A more detailed 
study of the parameter space reveals that the same is true for the right-handed down squarks. 

Light EW superpartners are a generic feature of the GGM parameter space. For example, 
we have seen that light Higgsinos in conjunction with light left or right-handed sleptons always 
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accompany the L and E boundaries of fig. 2 respectively. If fi > 0 we also expect a light wino 
throughout much of the parameter space, since the Higgs mass constraint selects out the 
neighborhood of M 2 ~ 0. This feature is especially important for low messenger scales, where 
the /i < 0 branch does not allow for light stops. Finally, due to the IR relation 

rn 2 A ~rnl 3 + n 2 , (5.1) 

the pseudoscalar may also be light. (Note that (5.1) holds anywhere in the parameter space, 
and will be strong test of GGM should these particles all be discovered.) The prevalence of all 
these light EW sparticles in GGM, often accompanied by decoupled colored sparticles, provides 
further motivation for dedicated Run II searches of direct EW superpartner production. 

Since the boundaries of the parameter space were determined by the combination of the 
Higgs mass, EWSB, and a tachyon, the GGM spectrum becomes especially predictive here. 
The tight connection between light EW states and the lightest possible mu 3 for a given M mess 
has a number of important consequences for collider searches. In particular, LEP bounds on 
EW states indirectly provide a lower bound on mu 3 . A future lepton collider such as ILC is 
expected to further probe a very large portion of the low mu 3 region of the GGM parameter 
space. 


5.3 Future directions 

We conclude by discussing some future directions. First, an important question is to what 
extent the constraints we derived here depend on the particular structure of the GGM boundary 
conditions at M mess . Here we briefly comment on more general scenarios: 


• The minimal extension of GGM relaxes the requirement of messenger parity, allowing 
for U (l)y H-tadpoles. This possibility was already discussed in [32]; see [49, 50] for 
explicit weakly coupled realizations. This breaks one of the sum rules in (1.1), leaving 
the residual ones: 

2 2 

m H d = ™ L 

rriq + 3 mfj — 9 m 2 D — 6 m 2 L + m 2 E = 0 (5.2) 

2 rriQ — 2>m\j + 3 m 2 L — 2m 2 E — 6 m 2 Hu = 0 

Consequently, an additional parameter must be added to the list (1.2). Full control over 
the resulting 8 dimensional parameter space might still be feasible by applying a similar 
strategy to the one we used here. The result can be interesting since the direct relation 
between mn u and rriL induced by (1.1) is now broken by the D-tadpole contributions. 
As a consequence the friction between large At and EWSB which was at the basis of our 
reasoning might be considerably alleviated and lighter stop masses could be viable. 
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• Adding a flavor blind mechanism to generate /r and B^ will generically break the sum- 
rules ( 1 . 1 ) in a model dependent way which cannot be parametrized by a reduced set of 
sum-rules like (5.2). Non-zero At may also be generated at M mess , which would obviously 
fundamentally alter the nature of the Higgs mass constraint. None of our conclusions 
can be then directly extrapolated to extended gauge mediation scenarios such as those 
in [12-28]. It may however be possible to perform a similar model independent analysis 
by making use of the framework developed in [34, 19]. 

• One can still focus on extended gauge mediation scenarios where A t is suppressed at 
M mess and is purely generated by RG evolution. This happens for example in solutions 
of the n/B^ problem which involve Higgs interactions with heavy singlets [34]. At can 
also be suppressed by an appropriate discrete R-symmetry [51]. In this context it would 
be interesting to account for the extra UV contributions to m 2 Hu and along the lines 
of what we have done here. As in the hypercharge D-term scenario, such contributions 
could alleviate the tension between EWSB and light stops, and possibly allow for the 
optimal point of m,h = 125 GeV with mg 3 ~ mjj 3 ~ \A t \/\/§ ~ 1 TeV. 

Another interesting direction for the future would be to study other aspects along the GGM 
parameter space which are not directly related to collider searches: 

• One of the peculiar features of GGM is that large At can only be achieved via large 
M 3 and therefore light stops require a careful tuning of the UV soft masses against the 
gluino RGE contribution. We expect this extra source of fine-tuning to play a significant 
role in the tuning measure. While we did not attempt to do so in this paper it would be 
interesting to quantify the tuning in GGM and comparing it against other UV complete 
gauge mediation models like the ones in [ 20 ]. 

• In order to compensate for the effect of heavy gluinos, the stops run tachyonic shortly 
above the weak scale [7]. Moreover, depending on the region of the parameter space, 
other UV tachyons are necessary in order to obtain EWSB in the IR. Thus the elec- 
troweak vacuum tends to be metastable. The estimates in [52] suggest that the vacuum 
decay is not a stringent constraint, but it would be interesting to perform a careful 
analysis in these GGM scenarios. 

• Cosmological bounds can also play an important role in GGM parameter space. In order 
to avoid gravitino overabundance and possible dangerous effects of the NLSP decays on 
the BBN products, inflation at particularly low temperature is required (see for example 
[31] and references therein for a discussion of the GGM cosmology). This bound on 
the reheating temperature can be evaded for example by adding tiny RPV couplings. 
However, it would be interesting to take it seriously and investigate in full generality the 
allowed cosmological scenarios in GGM. 
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Finally, let us mention that our procedure has some intrinsic limitations due to the theo¬ 
retical uncertainty in the Higgs mass computation. In this paper, we imposed rrih = 123 GeV 
in order to optimistically account for this uncertainty. It will be important to revisit this 
work after future improvements to the accuracy of the Higgs mass calculation, especially if 
these turn out to contribute negatively to rrih- (See for example the recent discussion in [53].) 
Aside from the usual corrections from higher orders and uncertainties in SM inputs like mt, 
the Higgs mass computation in GGM is particularly challenging due to the large hierarchies 
that are present in the colored spectrum. Perhaps the most acute example of such a situation 
is given by the Q boundary of our parameter space, where mjj 3> mg and also the gluino mass 
is very large. Using effective field theory techniques such as in [54], it would be interesting 
to have a better control on the Higgs mass computation in such a scenario where the lightest 
stop masses are realized in GGM. 
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A Validation plots 

In this section we will exhibit some plots validating the accuracy of the transfer matrix and 
the numerical algorithm for solving UV) = 0 and rrih = 123 GeV. 

Shown in fig. 16 are distributions of 5rriQ 3 and 5mu 3 across our entire scan of GGM pa¬ 
rameter space, which contains ~ 3 x 10 5 points. Here 5rriQ 3 and dmu 3 are the change in rag 3 
and mu 3 after running them up to the messenger scale using the transfer matrix and then 
back down to the weak scale using SoftSUSY. We see that the accuracy of the transfer matrix 
is very good, generally differing by less than ~ 50 GeV, and never differing by more than 
~ 200 GeV across the entire parameter space. These minor differences are due to effects not 
captured by the transfer matrix, primarily SoftSUSY’s iterative determination of Ms, and IR 
threshold corrections to 53 and yt . 

Shown in fig. 17 are rrih and \J -B^(UV) for every point in our GGM parameter scan. We 
see that the convergence on rrih is excellent, and the convergence on yj UV) is decent (99% 
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of the points have |y / .B /i (UV)| < 400 GeV). We also note in passing that the accuracy of 
the numerical scan is easily comparable to or larger than a naive estimate of the higher-loop 
B^(UV) expected from GGM. Given that the numerical scan was validated using the semi- 
analytic approach which assumed B^UV) = 0 , we do not expect that the small corrections to 
B^(UV) from GGM will make any difference to our conclusions. 




<5m Q3 (GeV) <5m U3 (GeV) 

Figure 16: Histograms showing the changes in toq 3 and mu 3 after running them up to the messenger 
scale using the transfer matrix and then back down to the weak scale using SoftSUSY, for every point 
in our GGM scan ( M mess = 10 7 , 10 11 , 10 15 GeV). 



yj (A4 mess ) (GeV) 



Figure 17: Histograms of rrih and \JB ll ( UV) across the entire GGM parameter scan. 


B Chargino/neutralino contribution to rrih 

In this appendix we will delve into the threshold correction to rrih from light winos and Higgsi- 
nos that greatly reduces the ^4-term required to achieve m.h = 123 GeV in the neighborhood 
of M -2 = 0 . 

Shown in fig. 18 is a plot of —At vs M 2 for a benchmark point extracted from our grid 
with M mess = 10 11 GeV and ( mQ 3 ,mu 3 ) = (1.5, 3.5) TeV. We see that A t varies by ~ 40% 
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for n > 0 (orange points), yet only varies by ~ 1.5% for /r < 0 (blue points). The variation for 
fj, > 0 comes in the form of a sharp decrease in the magnitude of At as we move from large M 2 
to small M 2 . This is characteristic of much of the parameter space, as we already saw in fig. 9. 

Clearly, the decrease in magnitude of At must be driven by a positive correction to the 
Higgs mass in the neighborhood of M 2 = 0, as shown in fig. 19. On the left, we have a plot 
of nth vs M 2 for the same benchmark of fig. 18. We see again the excellent convergence on 
nih = 123 GeV. Meanwhile in black we plot the same points, but hold At fixed at —2.8 TeV, 
which corresponds to the right-most point on the M 2 interval (indicated with a black star in 
fig. 18). We see that as we move towards M 2 = 0, if we don’t decrease the magnitude of A t to 
compensate for the effect at M 2 ~ 0, then the Higgs mass increases by as much as ~ 2.5 GeV. 

Alternatively one can see the same effect on the right plot of fig. 19. Here we show rrih vs 
At with all other soft parameters fixed to the values corresponding to the right-most point on 
the M 2 interval of fig. 18. We see that from A t = —2.8 TeV to A t = —2 TeV, rrih decreases by 
2.5 GeV. So we confirm that the origin of the variation in At is a ~ 2.5 GeV enhancement to 
rrih as we move to M 2 —» 0. 

It remains to isolate the origin of the 2.5 GeV threshold correction to rrih- Shown in fig. 20 
is a contour plot of rrih vs A 4 and M 2 , with all other parameters fixed to those of the right-most 
point on the M 2 interval. We see that the threshold correction is due to light charginos and 
neutralinos, and both light higgsinos and winos (i.e. both small /i and M 2 ) are required for the 
full effect. A similar effect was recently mentioned in the context of EFT calculations of the 
Higgs mass in [53]. We have further verified that this is the origin of the threshold correction, 
by direct calculation of the relevant one-loop diagrams as done in [46]. 


2.8 

^ 2.6 

> 

<u 

b 2.4 

b* 

1 2.2 
2.0 


Figure 18: Plot of —A t vs M 2 for the benchmark point indicated in the figure. Orange (blue) points 
correspond to /r > 0 (fj, < 0). The black star corresponds to rightmost point in the M 2 interval which 
has the biggest A t . 
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Figure 19: Left: Plot of rrih vs M 2 . Orange (blue) dots are grid points with \x > 0 (fi < 0). Black dots 
have At = —2.8 TeV as in the right-most point of the M 2 interval (the black star in fig. 18). Right: 
Plot of rrih vs —A t . We fix all the other parameters as in the right-most point of the M 2 interval in 
fig. 18. The dashed lines indicate the enhancement of the Higgs mass induced by the drop of M 2 which 
can be extracted again from fig. 18. 
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Figure 20: Contour plot of irih in the ( M 2 , fj) plane, with all other parameters fixed to those of the 
right-most point on the M 2 interval (the black star in fig. 18) and X t = A t — y t fi cot 13 kept fixed to —2.8 
TeV. Overlaid in orange (blue) are the actual values of M 2 and for the benchmark point in fig. 18 for 

/JL > 0 {/i < 0 ). 
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We conclude that the sharp decrease in the magnitude of At for /j > 0 is fully warranted. 
It corresponds to a few GeV threshold correction to the Higgs mass in the neighborhood of 
M 2 = 0 from light higgsinos and winos. It would be interesting to explore the implications of 
this threshold correction further in a more general context. Perhaps it could provide another 
motivation for EWKino searches at the LHC. 
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